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o. 

We study decoherence models for flavour oscillations in four-dimensional stochastically fluctuating 
space times and discuss briefly the sensitivity of current neutrino experiments to such models. 
We pay emphasis on demonstrating the model dependence of the associated decoherence-induced 
' damping coefficients in front of the oscillatory terms in the respective transition probabilities between 

flavours. Within the context of specific models of foam, involving point-like D-branes and leading 
to decoherence-induced damping which is inversely proportional to the neutrino energies, we also 
argue that future limits on the relevant decoherence parameters coming from TeV astrophysical 
neutrinos, to be observed in ICE-CUBE, are not far from theoretically expected values with Planck 
mass suppression. Ultra high energy neutrinos from Gamma Ray Bursts at cosmological distances 
can also exhibit in principle sensitivity to such effects. 

Oh 
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^ , I. INTRODUCTION 

The possibility that quantum gravity involve models with stochastic fluctuations of the associated metric field, 
' around some fixed background value, say flat Minkowski space time, may not be an unrealistic one. Such stochastic 
ON , models of gravity lead to observable consequences in principle, ranging from light cone fluctuations |l| to induced 
decoherence for matter propagating in such fluctuating space times 2]. Space time foam models naturally include 
such stochastic space time backgrounds. 

We describe here briefly the possible appearance of a stochastic space-time background by the D-particle foam 
model. It should be remarked that foam models do not consist exclusively of microscopic black holes, although 
this is the most popular model. The authors of [Sf] have provided another type of foam, inspired by brane-world 
scenarios. According to this model of D-particle foam, our brane world moves in a bulk space, punctured by point-like 
D-particles, which are point-like defects in membrane theory, characterised though by an infinity of super-Planckian 
string states. During such a motion, D-particles from the bulk cross our brane world, and interact with propagating 
matter in it, represented by open strings with their end-points attached on the brane. The interactions of string 
matter with the D-particle defects may not be universal among particle species. Due to electric charge conservation, 
electrically charged states, such as quarks or electrons may undergo only smooth ordinary scattering (recoil), while 
electrically neutral states, such as photons or neutrinos, may undergo topologically non trivial interactions in which a 
string state is split into two, with (some of) the corresponding end points are attached to the D-particle ("capture"). 
Capture includes recoil of the defect, as well as the induced back reaction onto space time, determined by means 
of conformal field theory methods of the associated non-critical string corresponding to the capture process 0. For 
heavy (Galilean) D-particles, the analysis leads to induced metric deformations of the (initially) Minkowski, say, flat 
background of the form: 

goo = —l,gij = Sij, gio = goi = -^ (1.1) 
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where i,j spatial (four-dimensional) indices on the brane world, with denoting a temporal direction, and Vi/c is the 
recoil 3-velocity of the D-particle defect, in units of the speed of light in vacuum, c, which by means of momentum 
conservation during the scattering process equals 

v, = g s Ah/M s . (1.2) 

In the last relation, Aki denotes the pertinent momentum transfer of matter, and M s /g s is the D-particle mass, with 
M s the string scale and g s < 1 the (weak) string coupling 0, 0] ■ We stress once more that the form (jl.ip is valid 
for Galilean (heavy) D-particles, relative to the momenta of the incident matter states. In general, for relativistic D- 
particles there are higher order corrections of Vi in the expression for the recoil-induced metric, leading to a covariant 
form in terms of the four-velocity = 7(1, v), with 7 the appropriate Lorentz factor [5[: 

9nu = f]nu + fi^u^Uy (1.3) 

where the scalar factor /($) depends on the details of the scalar dilaton $ configuration in the model. For instance, 
in the impulse approximation to D-particle recoil, the $ field is proportional to « fl X fl , as dictated by (logarithmic) 
conformal field theory considerations [5|, and /($) = 8 e (u At X /i ), with e — > + an appropriate regulator, linked to the 
world-sheet scale on account of closure of the conformal algebra Q. However, more general forms for the function 
/($) are possible, depending on the model. 

As discussed in Q , the above mentioned capture process may entail flavour oscillation as well as a non-trivial Fock 
vacuum structure of the vacuum experienced by the flavour states instead of the corresponding mass eigenstates. 
When considering statistical distributions of particles, one may have on average (ui) ~ 0, with only non trivial 
stochastic fluctuations (uiUj) = cr 2 8ij, a 1 7^ 0. The distribution of such D-particle velocities is a model dependent 
concept. In view of (|1.1|) . a stochastic statistical fluctuation of itj will also result in a stochastic fluctuation of the 
induced metric, and thus a stochastically fluctuating gravity theory. 

In we have studied the implications of such stochastic metrics in the simplified (but not unrealistic) case 
where the D-particle recoil velocities where along the direction of motion of the particle probe, thereby leading to 
an effectively two space-time dimensional problem, in which, however, the 7-matrix structures associated with the 
flavoured fermions (neutrinos) were kept four dimensional. For definiteness, in that work we have considered a model 
of Gaussian fluctuations, and we demonstrated that the associated oscillation probabilities of both fermion and boson 
flavoured particles exhibited exponential damping and decoherence, with the damping exponent scaling with the 
square of the time variable t. 

It is the purpose of this work to extend these results to more generic models of stochastically fluctuating metrics, 
living in a full-fledged four dimensional space time. In particular, we assume small random fluctuations , around 
the Minkowski background inverse metric r)^ , 

where for 7/ M „ we assume the sign convention (1, —1, — 1, —1) and W are defined as stochastic variables, obeying some 
statistical distribution, defining the model, such that (h^) — , but {h^hpu) are non trivial. 

The effect of neutrino oscillations will be studied in the above mentioned stochastic quantum gravity environment 
(|1.4p . Measurable quantities, such as transition amplitudes between neutrino flavours, will be obtained by averaging 
over the stochastic variables . As we will sec, the influence of the quantum gravity environment has a decoherence 
effect on the neutrino states, which is equivalent to a damping of the oscillations between neutrino flavours. 

The determination of the particular statistical distribution for the random fluctuations requires compete knowl- 
edge of the quantum gravity theory. In the absence of such a theory, the statistical distribution can not be fixed. 
Hence, our considerations below will be largely phenomenological, in an attempt to impose bounds on the relevant 
parameters by comparison to experiment. To this end, we mention that the allowed distributions must be such that 
the corresponding characteristic function yields finite, well-defined results for physical observables, such as the flavour 
oscillation probability. This criterion allows the use of distributions, such as the Cauchy-Lorentz (Breit-Wigner res- 
onance type), for which the mean (or odd moments) may not be well defined, while their variance (and higher even 
moments) diverge by themselves. For other distributions, such as the Gaussian of course the latter entities are 
well defined and small. It is these two kinds of distributions for the metric fluctuations that we shall concentrate 
upon in this work. Within the context of our D-particle foam model [3j, such distributions may characterise the 
ensembles of the D-particle recoil velocities. As we will see in the following sections, the choice of the distribution 
leads to different damping in the oscillation probability. For instance, the Cauchy-Lorentz case is characterised by 
exponential damping, with a linear power of time t, as in the Lindblad decoherence model @, 0], in contrast to the 
Gaussian model [2| , where the damping exponents are proportional to t 2 . 

In this work we shall be only concerned with stochastic gravity models, with fluctuations (|1.4[) about flat Minkowski 
space-time backgrounds. In such models we saw that the results for the flavour oscillation probabilities are identical 
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between fermions and bosons. The situation may be different, if the background space times are not Minkowski fiat. 
We shall come back to such more general situations in a forthcoming publication. Comparison of these two kinds, 
insofar as the order of magnitude of the expected effects, and their prospects of discovery in future experimental 
facilities, will be discussed in the concluding section of our article. 

The structure of the article is as follows: in section |TT] we set up the formalism for scalar flavoured particles, 
concentrating on the case of two generations for simplicity. We examine first the situation, in which only leading 
order (linear) terms in the deformation /i M „ are kept in the expressions for the pertinent oscillation probabilities, 
paying attention to demonstrate the dependence on the results, in particular the form of the damping exponents as 
functions of time, on the kind of the (statistical) distribution of the metric fluctuations. We also demonstrate that 
the inclusion of higher order corrections in hn V do not change the results qualitatively. In section IIIII we extend the 
analysis to two-generations of fermions, showing that the results concerning the oscillation probabilities and the effects 
of decoherence due to the stochastic metric fluctuations, are identical to the scalar case. Finally we comment briefly 
on a possible extension of the above results to the realistic three generation neutrino case, which however shall be 
considered in detail in a forthcoming publication. Conclusions and outlook are presented in section Hvl 



II. FLAVOUR OSCILLATIONS FOR SCALAR PARTICLES IN A STOCHASTIC METRIC 

BACKGROUND 



In this section we study massive scalar particles in the quantum gravity stochastic background (| 1 .4[) . The Klein- 
Gordon equation for a single scalar particle in a gravitational background, reads: 

(ff^VaV^ + to 2 ) = (2.1) 

where g a @ is the inverse metric tensor and V Q is a gravitational covariant derivative. 

A flat fluctuation h^" around the inverse Minkowski metric rf v is defined by Eq. 1|1.4[) . For simplicity, and motivated 
by the Galilean D-particle foam case (jl.ip , we may assumed that W v is independent of space-time coordinates [I?) • 
In such a case, the Christoffel symbols vanish, hence the Klein-Gordon equation in 3+1 dimensions reads: 

(g 00 d do + 2g 0i dodi + g'-'<),i), + to 2 ) <\> = (2.2) 

We seek for plane wave solutions of the form: 

fax,t) = 0(k,uj)e t{ujt - kx) (2.3) 

where k represents the momentum of the scalar particle in the Minkowski metric background. Note, that we have 
chosen one of the coordinate axes to lie along the direction of motion of the scalar particles, i.e. we have set = &3 = 
and ki = k. Plugging the expression (|2.3p into eq. (|2.ip leads to 

g 00 u; 2 ~2g 01 ku J + g 11 k 2 ~m 2 = (2.4) 

The two solutions uj(k), ui'(k) of the above equation are physically equivalent, since 

oj(-k) = -Lu'(k) . (2.5) 

Keeping the positive energy solution, we obtain the dispersion relation 

w = foo fc + ^oo V(9 01 ) 2 k 2 - g 00 ^ 1 ^ 2 - to 2 ) (2.6) 

We next proceed to discuss the flavour oscillations, by considering first the approximation of keeping at most linear 
terms in the stochastic fluctuations h^ v . 



A. Linear approximation for the stochastic metric fluctuations 



We assume for brevity two flavours for the scalar particles, written as a two-component scalar $ = (fa, fa)- The 
equation of motion is 



(g a0 y a y f 3 + m 2 ) $ = o 



(2.7) 
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where M 2 is the 2x2 positive definite (mass) 2 matrix, with eigenvalues m\ and m|. Note, that the flavour states are 
not eigenstates of M. For t = 0, we assume the production of a scalar particle of flavour a with density matrix: 

p(O) = |0 a )(^| (2.8) 

Taking into account that the energy eigenstates |/;) (I = 1,2) are related to the flavour states \<p a ) via the unitary 
transformation: 

|0 Q > = X^l/*) (^ a \=J2 U ^\ ( 2 - 9 ) 

i 3 

\.U) = E U UM (/il = E u n(<h\ ( 2 - 10 ) 
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we then obtain from Eqs. (|2.8[) and (I2.£)|) : 

pW^E 17 "^ ( 2 - n ) 

where we sum only over with a denoting the initial flavour state. The time evolution of the density matrix yields: 



P (t) =Y,U ai UZ i e i ^- u *\fi){fi\ (2.12) 



1,3 

Upon using Eq. (|2 . 10[) . one obtains the density matrix at time t, expressed in the basis of flavour states: 

P {t) = E U^U^U^e^-^p)^ . (2.13) 

h3,P,l 

The probability to find the scalar particle in a state of flavour f3 (/? ^ a) is 

Prob(a - /3) = E UaiU^U^-^ (2.14) 

i,3 

where the time dependent part is 

UaiU^U^Upie^-^* + U^U^U^Upie^-"^ (2.15) 
Since the perturbations ha V are stochastic variables, it is necessary to compute the following integral 

( e i(«i-a*)t) = J d 10 h F ( h ) e i(«i-a» 2 )* ( 2 16 ) 

where h = (h oa 7 h 01 , h 02 , h 23 , h 33 ), and .F(h) is a multi-variable probability density function. The measurable 
transition probability is found by averaging over the metrics 

(Prob(a - 0)) = E UaU^U^Uff^-"^) (2.17) 

In what follows we will assume that the stochastic variables are independent, or equivalently that: 

F(h) = f(h™)f(h nl )f(h° 2 ) x ... x f(h 23 )f(h 33 ) (2.18) 

where /(x) is a one-variable probability density function for the stochastic variable x. Additionally, we have assumed 
that the one- variable density function f(x) is the same for all W . We will approximate the energy difference lo\ — loi 
with a linear expansion of Eq. (|2.6p over the fluctuations h 

Wl (h) - w 2 (h) = a - ^ 00 - -/i 11 + 0(h 2 ) (2.19) 

where 



a = V fe2 + m i ~ V fc2 +m 2. fe= / = ( 2 - 20 ) 
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Note that in the above linear expansion only the elements h 00 and h 11 are involved. The integral (|2. 16[) then reduces 
to 

(e*" 1 " W2)t ) = J dh^dh 11 f(h ot) ) f(h lv ) exp {it (a - ah 00 /2 - bh u /2)} 

= $(-ai/2)$(-6t/2)exp(mt) (2.21) 

where 

/+oo 
e^ x f(x)dx (2.22) 
-oo 

is the characteristic function of the stochastic variable x with probability density function f(x). 

B. Gaussian fluctuations 

We will assume first a Gaussian distribution for the metric fluctuations, with density function 

2 / 2 

-x ja 

f(x) = -=, (2.23) 



where we have considered the case of zero mean value (fi = 0), and a is the corresponding standard deviation. The 
characteristic function for the Gaussian distribution is: 

i r + °° 

$(£) = -== / dx exp (i£x - x 2 /a 2 ) = exp (-£ V/4) (2.24) 

From Eq. ([2~2l]) we obtain 

(2.25) 



( e i(u™) t) = exp | ioi _ ( a 2 + 5 2) | 



Hence, the amplitude of the flavour oscillations has an exponential damping quadratic in time, that dependents on 
the particle's momentum. In the case of high energy particles (compared to their mass) one can make an asymptotic 
expansion for rrii/k << 1 

<e«— )*) . exp [iktA (l - } exp {-^A 2 } (2.26) 

where terms of order 0(rrii/k) 6 have been disregarded, and the parameter A is defined as 

A = ^ 2 « 1. (2.27) 

We remark at this point that self-consistency in the expansion in powers of rrii/k forces us to keep the forth-order 
correction term in the expression for the oscillation frequency in (|2.26p . since the damping is of order (rrn/k) A . The 
results are in full agreement with the corresponding effective two-dimensional case of 0] . A discussion on the prospects 
of observing such a decoherence effect in realistic physical systems of neutrinos, as well as on the dependence of the 
effect on the distance of the neutrino sources, especially in the case of astrophysical neutrinos, will be given in the 
concluding section of the article. 

C. Cauchy-Lorentz fluctuations 

We will consider next the case of Cauchy-Lorentz fluctuations 

fix) = -J—o (2-28) 
7r or + 7 Z 
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where 7 is the scale parameter, and the location parameter is zero. Note that the Cauchy-Lorentz distribution has 
no well defined mean, and in general odd higher-order moments, while its variance and higher-order even moments 
diverge. However, the characteristic function, defined as: 

i r°° 1 e^ x 

HO = - / dx \— 2 = exp (- 7 |£|) (2.29) 

yields finite results for the oscillation probabilities: 

(e *(^ 2 ) t) = exp | iat _ | ( | a | + (2 30) 

In this situation we obtain an exponential damping, with the exponent linear in time. An expansion in powers of 
rui/k gives 

(e^i-**)*) ~ exp UktA - 7*t|A|}. (2.31) 

where terms of order 0(rrii/k) have been disregarded. We did not write here the correction to the oscillation 
frequency, as this is of order (mj/fc) , whereas the damping is of order (mj/fc) 2 , and we are interested in the dominant 
effect only. The case is similar in form to the Lindblad decoherence [f| 0] • 



D. Beyond the linear approximation 

In this section we go beyond the linear approximation for Gaussian metric fluctuations, in order to check whether 

the previously derived damping of flavour oscillations in the linear-6, approximation is modified significantly. As we 
shall see, the main results remain valid. In this case we expand the energy difference oj\ — 0J2 up to quadratic terms 
in the fluctuations h 

wi(h) - uj 2 (h) = a + hd + h- Dh + C>(h 3 ) (2.32) 
with h = (h 00 , h 01 , h 11 ), d = (-a/2, 0, -6/2) and 

/3a/4 6/4 \ 

D = 6 (2.33) 
\ 6/4 -c/4 / 

The parameters a and 6 are given by Eq. (|2.20p . while c is defined as: 

fc 4 fc 4 

c = — o - — o ( 2 - 34 ) 

k l + mf k l + m| 

We next write the Gaussian probability density for the metric fluctuations in a compact form: 

F(h) = (2.35) 



where 



We wish to compute the integral 

e -h.B.h 



using the formula 



{e ^-^)t\ = / d 3 h , - e '("i-^)« (2.37) 



7T 3/2 /U-B _1 -U N 

d 3 h exp (-h • B • h + u • h) = -j— ^ exp f ) (2.38) 
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where the matrices B and u are defined as: 



B = H — itT> (2.39) 
u = itd (2.40) 



The matrix B can be written explicitly: 



^ - \ita ~\itb 
B = | 4j -itb j (2.41) 
-\itb " ±x + \itc 



Upon applying Eq. (|2.38p . we find: 



(e i( — 2) *> = fS|V /2 exp(-xW)exp(za<) 



V detB ) 
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where 

a 2 t 2 (4(a 2 + 6 2 ) - 4ia 2 t(b 2 - ac)) 



exp(-x(0) exp(iai) (2.42) 



X ™ 4a 4 i 2 (6 2 + Sac) - I6ia 2 t(3a - c) + 64 
P(f) = (l-ifecr 2 i)(16-4i(3a-c)cr 2 i+(6 2 + 3ac) ( T 4 t 2 ) (2.43) 



Expanding the exponent x(t) m powers of the small parameter a we obtain 

X (t) = ^-(a 2 + b 2 ) + ^(3a 3 + 2ab 2 -b 2 c) 
16 64 

- — (9a 4 + 7a 2 6 2 + b 4 - 2ab 2 c + b 2 c 2 ) + 0(a 8 ) (2.44) 
256 

The leading term of the above expansion corresponds to the main effect which is the exponential damping of particle 
oscillations. Note that it is identical with the one found in the linear approximation above. The next to leading order 
term is purely imaginary and modifies slightly the oscillation term e mt . The factor P(t) has a subleading contribution 
of the form: 

\P(t)\ = 16 + ^(9a 2 + 18b 2 + c 2 ) + 0(a 6 ) (2.45) 

The reader should compare the results with the effectively two-dimensional model of 0). The results are similar, 
independently of the representation of gravity fluctuations and the number of spatial dimensions. 



III. FERMIONS IN A STOCHASTIC METRIC BACKGROUND 



In this section we would like to extend the above results to incorporate particles with spin (fermions), which from 
a phenomenological point of view is more interesting, in view of the potential application to the physics of neutrino 
oscillations. However, as we shall demonstrate below, the results will be similar to the bosonic case, as far as the main 
features of decoherence damping is concerned. This is attributed to expanding about a flat Minkowski background 

C3J). 



A. Dirac fermion dispersion relation 

In this section, we will consider Dirac fermions, but the results presented are the same for Majorana fermions, as the 
only difference is that the latter are described by self-conjugate fields, which does not affect the dispersion relation. 
We remark that this difference would play a role only if we considered the MSW effect Q, for the propagation of 
neutrinos in matter media, as in that case only one Weyl spinor component of the fermion would be affected. We shall 
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consider this case, along with more general situations involving expansions about non flat backgrounds in a future 
publication. 

We review here the basic elements we will need to describe fermions in a curved background. We shall follow the 
formalism of 0] and in references therein, where we refer the reader for details. At each point M of space time, the 
vierbeins e a = <9 Q M span the flat tangent space time, and are related to the inverse metric by 

<e v f>rf*=gr (3.1) 

The gamma matrices 7° which follow are defined on the flat tangent space time and satisfy {7", 7 /3 } = 2?7 Q/3 , and we 

also define ct q/3 = i[j a , j }/2. 

The Dirac equation in a curved background is 

{ij a V a - m) V = (3.2) 

where 



(3.3) 



and is the covariant derivative. In the present situation, we are interested in a constant and homogeneous stochastic 
metric, such that the vierbeins do not depend on space time coordinates, and the Christoffcl symbols vanish. As a 
consequence, the Dirac equation reads 

(ij a e^ - m) i> = (3.4) 



In order to find the dispersion relation for the fcrmion, we multiply the Dirac equation (13.41) by the complex conjugate 
operator (— ij^e^d v — m) to obtain 

= Q{7 a ,7^K e^d v +rrAi, 

= (g^d^ + m 2 )^ (3.5) 

which is similar to the equation (|2.ip for a boson, since covariant derivatives are simple derivatives in our case. As a 
consequence, the spin does not play a role in the situation where the background metric is flat, and the previous results 
derived for a boson apply to fermions. It appears that damping is a general result of quantum gravity fluctuations. 



B. Two- flavour fermion oscillations 



We consider two Dirac fermions ipe,^^, written as an eight-component fermion 'J, which are coupled by a Dirac 
mixing mass matrix, and are described by the equation of motion 

(ij a V a - M) * = 0, (3.6) 

where the mass matrix in flavour space reads: 



M = e ^ (3.7) 
1 m Rl , m u 1 



In order to involve the mass eigenstates ipij^, one performs the following rotation in flavour space 

(3.8) 



where 



ip e \ _ / cos 8 sin 9 \ I ipx 
iftjj, J ~ \ — sin 8 cos 9 J yip2 

tan(20) = 2WeM (3.9) 
- m e 

As explained in (To| , the Lagrangian describing the fermions tpi and ^2 is then the sum of two independent Lagrangians, 
such that the equations of motion for ipi and ip2 are 

(?7 Q 2? Q - mj) ^ - i = 1,2, (3.10) 
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where mi + m 2 = m e + and m\mi = nriera^ — m\^. Hence, this leads to the study of two individual fermions. 
The results are then similar to those derived for bosons, as the fermion dispersion relation is identical with that of 
bosons, as can be readily seen from Eq. (|3.5p . 

Thus, using Eq. (|2.14[) and Eqs. (|2.25p , (|2.31[) and (|3.8[) . one obtains the formula for the transition probability 
between flavours, in the linear-order approximation for the /i-metric fluctuations: 

(Prob(a -» p)) = i sin 2 (26>) (l - e" x(t) cos(at)) , (3.11) 

where a ^ (3 and 

• x{t) = ^rtr ( fl2 + f° r the Gaussian case; 

• x{t) = ~2~(\ a \ + H) f° r the Cauchy-Lorentz case. 

We note that, in the case where there is no metric fluctuation (for %(i) = 0), the probability (|3.11j) leads to the known 
result | sin 2 (2#)[l — cos(af)] [ll|. Taking into account metric fluctuations, we see that the limit t — > oo leads to the 
stationary probability h sin 2 (29), which is independent of the fermion energy. 

C. Three-flavour fermion oscillations 

It is not difficult to generalize our results in the case of three flavours. As we will see the exponential damping 
remains. 

The corresponding 3x3 Dirac mass matrix M, defined in analogy with Eq. (|3.7[) . can be diagonalized with a 
unitary transformation U a i. 

\il> a ) = X) Uai\il>i) (a = e,n,r) (3.12) 

i 

and has three positive eigenvalues m, (i — 1,2,3). We note for completeness that for the U matrix we use the 
parametrization of (llj . 

The probability to find a fermion, with initial state a, in a final state with flavour (3 ^ a is given by Eq. (|2.17p : 

(Prob(a -f 13)) = U ai UZ j U$ i Ufi j (e« u *- u #). (3.13) 

As in the two-flavour case, the probability (|3.13[) has a time-independent part, which is 

I U al | 2 | Up | 2 + | U a2 | 2 | Up 2 | 2 + I U a3 | 2 | U 03 | 2 , (3.14) 
whereas its time-dependend part is 

UaiU^U^U^ie^-^) + U^U^UfcUpsie^-"^) + U^U^U^U^e^-^) + cc. (3.15) 
We then average over the stochastic fluctuations , and calculate the three elements 

/gi(u)i— tt) 2 )t\ /gi(ct>2— w 3 )t\ /gi(tt) 3 — wi)t\ 

Using our previous results (see Eqs. (|2.26[) and (|2.3ip ), we obtain for Gaussian fluctuations 

* «p{«W« (l - j exp {-^2!a» } , (3-16) 

and for Cauchy-Lorentz fluctuations 

(e^*-^)*) ~ exp {tfctAy - 7 fci|Ay |}, (3.17) 

where 

2 2 

mf - rrij 



= ^2fc^ « L ( 3 - 18 ) 
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The reader is reminded at this point that, as with (|2.26p . consistency in the expansion in powers of rrii/k forces us 
to keep the forth-order correction in rrii/k in the oscillation frequency (|3.16[) of the Gaussian model, given that the 
damping is of order (rni/fc) 4 . On the other hand, in the Cauchy-Lorentz case (|3.17p we disregarded such a forth order 
correction, since the corresponding damping is of order {rrii/k) 2 . 

We now notice that the l/i?-dependent exponential damping in the case (|3.17[) may be translated [HI, [lij] as 
implying a finite life time Tub for the probe in the laboratory frame: 

where r Vi is the life-time of the neutrino probe in the rest frame of the massive neutrino. It is interesting therefore to 
use our microscopic models in order to identify a possible quantum-gravitational origin of a finite lifetime of neutrinos. 
We shall come back briefly to this point later on the article, when we discuss lower bounds on the decoherence-induced 
neutrino lifetime (|3.19p from (current and future) experiments. 



IV. CONCLUSIONS AND OUTLOOK 



In this concluding section we would like to discuss briefly the above results in light of the prospects for observing 
the above effects in current or upcoming neutrino facilities. A more detailed comparison with experimental data and 
derivation of bounds will appear in a forthcoming publication. The main results of our work was the exponential 
decoherence-induced damping in the oscillation probability, which for the case of Gaussian fluctuations, with variance 
a 2 , is given by (|3.16p , (|3.18p . while for the case of Cauchy-Lorentz distribution by (|3.17p . 

For realistic models of quantum gravity, we first observe that in the Gaussian model, the damping exponent is 
much more suppressed than the corresponding one in the Cauchy-Lorentz (Lindblad-like-time-scaling) case. This is 
due to the extra suppression factor (m\ — m\)/E appearing in the exponent of the Gaussian model. In order to 
make direct comparison with experiment, it is essential to use as concrete quantum gravity models as possible. To 
this end, the D-particle foam model [H turns out to be very useful. According to this model, the stochastic metric 
fluctuations are induced, as explained in the introduction, by means of the recoil of the D-particle defect during its 
topologically non-trivial interaction/capture process with matter. Considering a gas of D-particles, it is natural to 
consider a Gaussian model for the distribution of the respective (spatial) recoil velocities, Ui, which in turn induce 
the metric fluctuations (jl.ip . for heavy defects. In such a case, the distribution variables are the dimensionless ratios 
Ui/c, which however are restricted to be less than one in magnitude, as are not allowed to exceed the speed of light 
in vacuum c, as requested by string theory. In such a case, the relevant integrals in the previous sections, e.g. (I2.24p . 
have to be understood to be restricted to the range [—1,1]. 

In this sense, the above results can only be viewed as an idealisation of the situation characterising the D-foam 
model, which however, is a pretty good approximation of reality if the variance of the Gaussian distribution is pretty 
small, which is the case in the D-foam model. Indeed, taking into account that a 2 in such a case might naturally 
be expected to have an order of magnitude dictated by the square of (|1.2p . since after all a 2 describes the variance 
of D-particle recoil velocities, which was assumed small, and the typical order of the latter is given by (|1.2p . If one 
assumes that the momentum transfer of a particle probe, during its capture by a recoiling D-particle spacetime defect, 
is of the order of the particles momentum itself, then we may have for the variance a 2 the estimate: 

2 _ /n I a „2 



This yields a damping exponent (|3.16p of order (in units h = c = 1): 

exo(-ft 2 t 2 ) n 2 - gMl^ f42 ) 

CX P \ "Gauss r ) ) "Gauss— ^ ]\,f2 ' * • ' 

For the case of Cauchy Lorentz distribution of D-particle velocities (adopted appropriately, as in the Gaussian 
case, to incorporate velocity variables Ui which do not exceed the speed of light in vacuum c), one observes that the 
parameter 7 provides a characteristic scale for D-particle velocities, since, if they exceed that scale the distribution is 
diminished significantly. On account of (|1.2p . this leads to the assumption that a natural order of magnitude for 7, in 
the context of the D-particle foam model, is 7 ~ g s E/M s , with E the energy of the particle probe. The corresponding 
decoherence-induced damping exponent has then the form: 

ex P (-n CL t) , n CL = ^' 2M 3 ' . (4.3) 
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Searches for Lindblad decoherence [12| , using the latest neutrino experimental data, have bounded the respective 
coefficients in a stringent way. For two-flavour models, the parametrization used in [12j for the decoherence-induced 
Lindblad- type damping coefficients is: 



exp (-7 t) 



* = 70 ( Sv) 



(4.4) 



with the following bounds provided by means of combining atmospheric, solar-neutrino and KamLand data [lz 



70 < 0.67 x 10~ 24 GeV 
7o < 0.47 x 10~ 20 GeV 
7o < 0.78 x 10~ 26 GeV 



n = 
n = 2 
n = — 1 



(4.5) 



It should be remarked that all these bounds should be taken with a grain of salt, since there is no guarantee that in a 
theory of quantum gravity 70 should be the same in all channels, or that the functional dependence of the decoherence 
coefficients 7 on the probe's energy E follows a simple power law. Complicated functional dependencies j(E) might 
be present in general. 

To compare with our model above (|4.3I) . we observe that the resulting damping coefficient is independent of the 
probe's energy E, as a result of the -E-dependence of the coefficient 7. For such constant decoherence coefficients, the 
analysis of 12j, (|4.5[) . yields the following bound on M s /g s : 



M, 



> 0.74 x 10 24 max, 



GeV' 



GeV 



(4.6) 



where max^ j indicates the maximum mass difference among neutrino flavours, and we assumed that the quantum 
gravi ty p arameters are the same for all flavours, which is certainly the case of the D-particle foam model. Recent 
data indicate that Am? G (10~ 23 - 10~ 21 ) GeV 2 , from which gl| implies: 



M, 
9s 



> 740 GeV 



(4.7) 



The expected minimal value of the string mass scale M s is, of course at least a few TeV, for which the string 
coupling must be very weak in order to provide realistic string phenomenology. For couplings of order g s < 1/2 the 
phcnomenologically correct scale is close to four-dimensional Planck scale, M s ~ 10 18 GeV. The above considerations, 
therefore, imply that the currently available neutrino data do not have the sensitivity to probe realistic D-particle 
foam models (the Gaussian Models decoherence is much more suppressed than the Cauchy-Lorentz one, as already 
mentioned) . 

One may reverse the logic, and try to consider bounds in as much model independent way as possible, using the 
results (|3.16p . (|3.17[) . without reference to any explicit model for the parameters a 2 and 7. In such a case, one may 
assume that these parameters are probe-energy independent, and try to bound their values, by comparing with data. 
We postpone such a complete analysis for a forthcoming work. Here, however, we note in brief that in the Cauchy- 
Lorentz case with constant scale-parameter 7, the resulting decoherence coefficient corresponds to the l/i?-dependent 
case, n = — 1 in (|4.4[) . for which, on account of (|4.5p . one obtains the following bound on 7: 



7 < 10" 



(4.8) 



on account of the above-mentioned measured neutrino mass differences. 

A final comment, concerns the order of the damping exponents in (|3 . and (|3.17[) for astrophysical neutrinos. To 
answer such a question it is imperative to know the energy dependence of the respective damping coefficient. Since 
the decoherence effects depend on the actual time of flight of neutrino t, the effects are maximised for extraterrestrial 
neutrinos, coming from extragalactic sources. Potential limits on Lindblad decoherence using high e nerg y (TeV scale) 
astrophysical (anti)neutrinos from ICE-CUBE in the future have been analysed in [l4[ and also in jl5l | In the ICE- 
CUBE case 14]. the inverse energy n = — 1 decoherence coefficient 70 in the terminology (|4.4[) . is found to be (at a 90% 
CL) 70 < 10 -34 GeV, which will improve the existing sensitivity by eight orders of magnitude, implying, for a constant 
Cauchy-Lorentz distribution, the limit 7 < 10 -13 . Such sensitivities are not far from naturally expected values of 7 in 
space-time foam models involving heavy (Planck-mass M s /g s ~ 10 19 GeV) D-particles and TeV-momentum transfers 
during the capture of TeV- en ergy (anti)neutrino matter by the recoiling D-particle. Moreover, for completeness, we 
mention that, as discussed in the 90% CL bound on the inverse-energy decoherence from ICE-CUBE will imply, 



12 



according to (|3 . 19[) the existing bounds on the electron-antineutrino life time Tv ts /mv e > 10 34 GeV -2 , improving by 
four orders of magnitude the existing bounds from solar neutrinos. 

Finally, we mention that in [n| it was argued that ultra-high energy neutrinos, with energies 10 17 — 10 19 eV can be 
emitted by Gamma Ray Bursts (GRB), which actually carry a significant fraction of the GRB energy. Most of GRBs 
lie at cosmological distances corresponding to distances z > 1, i.e. larger than 10 27 m. To study self-consistently such 
cases, one needs of course to consider the propagation of the neutrino in a Robertson- Walker background, about which 
one could expand the space-time metric fluctuations. We hope to come to a study of such issues in a future work. 
However, to get a preliminary idea it suffices to ignore the expansion of the Universe and consider the Cauchy-Lorentz 
decoherence distribution as a pilot case. For such distances and energies, the relevant exponent (|3. 17|) of a Cauchy- 
Lorentz decoherence model with constant (probe-energy-independent) 7 < 10~ 5 (c.f. (|4.8p ). becomes roughly of order: 
5 x 10~ 37 x t, where we took into account that the maximal of the neutrino mass difference is max; .,- Am^- ~ 10 -3 . 
Such damping becomes of order one for distances of order L ~ 10 27 m, i.e. at cosmological distances of GRBs. Of 
course, the expansion of the Universe, will modify such results, nevertheless what this preliminary exercise showed is 
that high energy astrophysical neutrinos can indeed constitute sensitive probes of decoherence models [3, EH • We 
hope to be able to come back to a detailed discussion of such issues in the future. 
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